Effective Dissipation and Turbulence in Spectrally Truncated Euler Flows 
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ABSTRACT: A new transient regime in the relaxation towards absolute equilibrium of the conser- 
vative and time-reversible 3-D Euler equation with high-wavenumber spectral truncation is charac- 
terized. Large-scale dissipative effects, caused by the thermalized modes that spontaneously appear 
between a transition wavenumber and the maximum wavenumber, are calculated using fluctuation 
dissipation relations. The large-scale dynamics is found to be similar to that of high-Reynolds 
number Navier-Stokes equations and thus to obey (at least approximately) Kolmogorov scaling. 

PACS numbers: 47.27.Eq,05.20.Jj, 83.60.Df 



Turbulence has been observed in inviscid and conser- 
vative systems, in the context of (compressible) low- 
temperature superfluid turbulence 0, 0, 0- This be- 
havior has also been reproduced using simple (incom- 
pressible) Biot-Savart vortex methods, which amount to 
Eulerian dynamics with ad hoc vortex reconnection 
The purpose of the present letter is to study the dy- 
namics of spectrally truncated 3-D incompressible Eu- 
ler flows. Our main result is that the inviscid and con- 
servative Euler equation, with a high-wavenumber spec- 
tral truncation, has long-lasting transients which behave 
just as those of the dissipative (with generalized dissi- 
pation) Navier-Stokes equation. This is so because the 
thermalized modes between some transition wavenumber 
and the maximum wavenumber can act as a fictitious 
microworld providing an effective viscosity to the modes 
with wavenumbers below the transition wavenumber. 

We thus study general solutions to the finite system of 
ordinary differential equations for the complex variables 
v(k) (k is a 3 D vector of relative integers (^1,^2,^3) 
satisfying sup Q \k Q \ < fc max ) 

d t v a (k,t) = --V a 0j(k)^2vf3(p,t)v 7 (k- p,t) (1) 

p 

where V a ^ — kpP ai + kjP a p with P a p = S a p — k a kj3 /k 2 
and the convolution in Q is truncated to sup Q \k a \ < 

"-max j 

sup Q \p a \ < fcmax and sup Q \k a — p a \ < k 

max ■ 

This system is time-reversible and exactly conserves 
the kinetic energy E = E{k,t), where the energy 
spectrum E(k, t) is defined by averaging v(k', t) on spher- 
ical shells of width Ak — 1, 

E(k,t) = \ E Iv(k',i)| 2 . (2) 

fc-Afc/2<|k'|<fc+A/c/2 

The discrete equations Q are classically obtained 
by performing a Galerkin truncation (v(k) = for 
sup a |fc Q | < fcmax) on the Fourier transform v(x, t) — 



J^v(k, i)e lkx of a spatially periodic velocity field obey- 
ing the (unit density) three-dimensional incompressible 
Euler equations, 

d t v + (v • V)v = -Vp , 

V • v = . (3) 

The short-time, spectrally-converged truncated Eulerian 
dynamics Q has been studied ja, Q to obtain numer- 
ical evidence for or against blowup of the original (un- 
truncated) Euler equations (0- We will study here the 
behavior of solutions of £Q) when spectral convergence 
to solutions of © is lost. Long-time truncated Eulerian 
dynamics is relevant to the limitations of standard simu- 
lations of high Reynolds number (small viscosity) turbu- 
lence which are performed using Galerkin truncations of 
the Navier-Stokes equation || . 

Equations Q are solved numerically using standard 
pseudo-spectral methods with resolution N. The so- 
lutions are dealiased by spectrally truncating the modes 
for which at least one wave-vector component exceeds 
N/3 (thus a 1600 3 run is truncated at fc ma x = 534). 
This method allows the exact evaluation of the Galerkin 
convolution in (JIJ in only N 3 \ogN operations. Time 
marching is done with a second-order leapfrog finite- 
difference scheme, even and odd time-steps are period- 
ically re-coupled using fourth-order Runge-Kutta. 

To study the dynamics of Q, we use the so-called 
Taylor-Green ^(J single-mode initial condition of © 
u TG — sin x cosy cos z , v TG = — u TG (y, —x, z), w TG = 0. 
Symmetries are employed in a standard way to re- 
duce memory storage and speed up computations. Runs 
were made with N = 256, 512, 1024 and 1600. 

Figure n displays the time evolution (top) and resolu- 
tion dependence (bottom) of the energy spectra. Each 
energy spectrum E(k,t) admits a minimum at fc = 
fcth(i) < fcmax, in sharp contrast with the short-time 
(t < 4) spectrally converged Eulerian dynamics (data not 
shown, see pHllj). For fc > fcth(i) the energy spectrum 
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FIG. 1: Energy spectra, top: resolution 1600 3 at t = 
(6.5,8,10,14) (o, +, o, *); bottom: resolutions 256 3 (circle 
o), 512 3 (triangle A), 1024 3 (cross x) and 1600 3 (cross +) at 
t = 8. The dashed lines indicate k~ 5 ^ 3 and k 2 scalings. 



obeys the scaling law E(k,t) = c(t)k 2 (see the dashed 
line at the bottom of the figure). The dynamics thus 
spontaneously generates a scale separation at wavenum- 
ber fcth(£)- Figure also shows that fcth slowly decreases 
with time. For fixed fc inside the k 2 scaling zone E(k, t) 
increases with time but E(k, t) decreases with time for k 
close (but inferior) to fcth(i)- 

The traditionally expected [Gl.ll2j asymptotic dynamics 
of the system is to reach an absolute equilibrium, which is 
a statistically stationary exact solution of the truncated 
Euler equations, with energy spectrum E(k) = ck 2 . Our 
new results (see figure \Q show that a time-dependent 
statistical equilibrium appears long before the system 
reaches its stationary state. Indeed, the early appearance 
of a k 2 zone is the key factor in the relaxation of the sys- 
tem towards the absolute equilibrium: as time increases, 
more and more modes gather into a time-dependent sta- 
tistical equilibrium which itself tends towards an absolute 
equilibrium. 

Since the total energy E is constant, the energy dissi- 
pated from large scales into the time dependent statistical 
equilibrium is given by 

E th (t) = J2 ^(M) ■ (4) 
fc th (t)<fc 



The time evolutions of fcth and Eth are presented on figure 
12 The figure clearly displays the long transient during 
which, for all resolutions, fcth decreases and E t h increases 
with time. Note that, at all times, fcth increases and Eth 
decreases with the resolution. 

Since the energy of the time-dependent equilibrium 
increases with time, the modes outside the equilibrium 
lose energy. The presence of a time-dependent equilib- 
rium thus induces an effective dissipation on the lower fc 
modes. 

We now estimate the characteristic time of effective 
dissipation T(fcd) of modes fcd close to fcth(^) by assuming 
time-scale separation and studying, at each time t, the 
relaxation towards the time-independent absolute equilib- 
rium characterized by Eth{t) and fc max . The existence of 
a fluctuation dissipation theorem (FDT) [l^L fl3 | ensures 
than dissipation around the equilibrium has the same 
characteristic time-scale as the equilibrium correlation 
functions (£> Q (k, t)-0a(k', 0)) (brackets denote equilibrium 
statistical averaging over initial conditions £>^(k',0)). 
Defining this time scale tq as the parabolic decorrela- 
tion time 

7ga tt <fi a (k,i)^(k / ,0))| t=0 = <i5a(k,0)t) i 8(k' > 0)) J (5) 

time translation invariance allows to express the second 
order time derivative as — (d t v a (k, t)d t 'Vp(k! , t')) i t=t , =0 • 
Using expression Q for the time derivatives reduces 
the evaluation of tq to that of an equal-time fourth- 
order moment of a gaussian field with correlation 
< v a (k,t)vp(-k,t) >= AP a0 (k) where A = 
-Eth/(2fc m ax) 3 ■ The only non- vanishing contribution is 
a one loop graph 0, 0] . The correlation time rc associ- 
ated to wavenumber fc is found in this way [l^j to obey 
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the simple scaling law 



C 



TG 



(6) 



where C = 1.43382 is a constant of order unity. The 
time-scale tq is the eddy turnover time at wavenumbcr 
fcth- Because of Kolmogorov (K41) behavior (see be- 
low) the evolution of E t h is governed by the large-eddy 
turnover time. The assumption of time-scale separation 
made above is thus consistent. 




FIG. 3: Time evolution of the ratio fcth/fcd at resolutions 256 3 
(circle o), 512 3 (triangle A), 1024 3 (cross x) and 1600 3 (cross 

+) • 



This strongly suggests to introduce an effective gen- 
eralized Navier-Stokes model for the dissipative dynam- 
ics of modes k close to feth(t). To wit, we make the 
Ansatz e(k,t) = D\k\E(k,t), where v = ^/E^/C and 
e(k,t) = —dE(k,t)/dt is the spectral density of energy 
dissipation 



e(t) 



cLEthft) 
At 



Assuming that this dissipation takes place in a range of 
width akd around fed , we estimate the total dissipa- 
tion s ~ Pk ( iE(k ( {)ak ( i. This, together with -E(fcd) ~ 
IX yields the relation 



klE^/kl 



1/4 



E\ 



3/2 
th , 



1.3/4 

m -V" 



(8) 



The consistency of this estimation of effective dissipation 
with the results displayed in figure |3 requires that fed ~ 
/c t h- The ratio fcth/^d is displayed on figure [3J ft is seen 



to be of order unity and is reasonably constant in time 
and resolution independent (at least for TV > 256). 

Thus the small-scale modes between fc t h and fc max act 
as a fictitious thermostat providing, via the FDT, an ef- 
fective viscosity to the large-scale modes with wavenum- 
bers below fcth- Note that spontaneous equilibration hap- 
pening in conservative isolated systems, such as the one 
studied in the present letter, should not be confused with 
equilibration resulting from interaction with the ther- 
malized degrees of freedom of the molecules constitut- 
ing a physical fluid. Indeed the reversible dynamics of 
the isolated system Q spontaneously generates both the 
wavenumber at which the fictitious thermostat begins 
and its temperature. 




FIG. 4: Temporal evolution of, top: energy dissipation e (left 
vertical axis) and e~ 1,/3 (right vertical axis); bottom: k~ n 
inertial range exponent n at resolutions 256 (fit interval 2 < 
k < 12, circle o), 512 3 , (fit interval 2 < k < 14, triangle A), 
1024 3 (fit interval 2 < k < 16, cross x and 1600 3 (fit interval 
2 < k < 20, cross +). 



The previous results indicate scale separation between 
conservative large-scale and dissipative small-scale dy- 
namics. Furthermore the scale separation increases with 
resolution. This strongly suggests that large-scale be- 
havior may be identical to that of high- Reynolds number 
standard Navier-Stokes equations, which is known |8( to 
obey (at least approximately) K41 scaling. 

The energy dissipation rate Q shown on figure0|(top, 
left axis) is in good agreement with the corresponding 
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data for the Navier-Stokes TG flow (see reference 
figure 7 and reference figure 5.12). Both the time for 
maximum energy dissipation t max ~ 8 and the value of 
the dissipation rate at that time e(t max ) ~ 1.5 are 
in quantitative agreement. Furthermore the long-time 
quasi-linear behavior of e -1 / 3 (shown on right axis) is 
compatible with K41 self-similar decay s(t) ~ L\t~ z . 

A confirmation for K41 behavior around i ma x is dis- 
played on figure 0] (bottom) . The value of the inertial- 
range exponent n, obtained by low-fc least square fits 
of the logarithm of the energy spectrum to the function 
c—n log(fc), is close to 5/3 (horizontal dashed line) when 
t ~ tmax- The —5/3 exponent is also shown as the left 
dashed line on bottom of figure where the dissipative 
effects can be traced back to the energy spectrum de- 
creasing faster than fc -5 / 3 at intermediate wavenumbers. 

The mixed K41/absolute equilibrium spectra have al- 
ready been discussed in the wave turbulence literature 
(e.g..|lfij) and have more recently been studied in connec- 
tion with the Leith model of hydrodynamic turbulence 
[l7j . In this context, small-scale thermalization may have 
some bearing on the so-called bottleneck problem if the 
dissipation wavenumber approaches fc max . 

Note that the dynamics of spectrally truncated time- 
reversible nonlinear equations has also been inve stig ated 
in the special cases of 1-D Burgers-Hopf models [l8j and 
2-D quasi-geostrophic flows |l9j . A central point in these 
studies was the nature of the statistical equilibrium that 
is achieved at large times. Several equilibria are a priori 
possible because both (truncated) TD Burgers-Hopf and 

2- D quasi-geostrophic flow models admit, besides the en- 
ergy, a number of additional conserved quantities. The 

3- D Euler case is of a different nature because (except for 
helicity that identically vanishes for the flows considered 
here) there is no known additional conserved quantity |£| 
and the equilibrium is thus unique. The central prob- 
lem in truncated 3-D Eulerian dynamics is therefore the 
mechanism of relaxation towards equilibrium, as studied 
in this letter. 

In summary, our main result is that the spectrally trun- 
cated Euler equation has long-lasting transients behaving 
just like those of the dissipative Navier-Stokes equation. 
The small-scale thermalized modes act as a fictitious mi- 
croworld providing an effective viscosity to the large-scale 
modes. These dissipative effects were estimated using a 
new exact result based on Fluctuation Dissipation rela- 
tions. Furthermore, the solutions of the truncated Euler 
equations were shown to obey, at least approximately, 
K41 scaling. In this context, the spectrally truncated 
Euler equations appears as a minimal model of turbu- 
lence. 
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